Lifetime statistics in transitional pipe flow 
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Several experimental and numerical studies have shown that turbulent motions in circular pipe 
flow near transitional Reynolds numbers may not persist forever, but may decay. We study the 
properties of these decaying states within direct numerical simulations for Reynolds numbers up to 
2200 and in pipes with lengths equal to 5, 9 and 15 times the diameter. We show that the choice 
of the ensemble of initial conditions affects the short time parts of lifetime distributions, but does 
not change the characteristic decay rate for long times. Comparing lifetimes for pipes of different 
length we notice a linear increase in the characteristic lifetime with length, which reproduces the 
experimental results when extrapolated to 30 diameters, the length of an equilibrium turbulent puff 
at these Reynolds numbers. 
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I. INTRODUCTION 



In several shear flows such as plane Poiseuille, plane 
Couette [Tl, ^ and also pipe flow % i, i, i, Q turbu- 
lent dynamics is observed for flow speeds where the lam- 
inar profile is still stable against infinitesimal perturba- 
tions. In such a situation a finite amplitude perturbation 
is required to drive the system from laminar to turbulent 
flow [3, Q , and one might expect that also for the con- 
verse process, returning from the turbulent dynamics to 
the laminar one, a sufficiently large perturbation on top 
of the turbulent dynamics should be required. Several 
observations in direct numerical simulations and exper- 
iments show, however, that turbulent motion returns to 
the laminar flow suddenly and without any noticeable 
precursor or perturbation ^9, 10, 11, 12,]. From the point 
of view of nonlinear dynamics such a behavior suggests 
that the turbulent state does not correspond to a closed 
off turbulent attractor but rather to an open turbulent 
chaotic saddle [§, [ll|, [11] . One can then assign to each 
initial ffow state a lifetime, i.e. the time it takes for this 
state to return to the laminar profile. The lifetime is 
a valuable observable that has previously been used to 
extract information about states on the border between 
laminar flow and turbulence [l3, [3 ■ ^'^^^ hcic use 
it to extract information about the turbulent dynamics 
itself, thereby extending the work reported in Ref. [ll[. 

Experiment and simulations show that neighboring 
trajectories can have vastly different lifetimes, so that the 
lifetime is rather unpredictable and depends sensitively 
on the initial perturbation, see e.g. [lO, [HI, El]- This 
strong sensitivity on initial conditions is consistent with 
observations on other transiently chaotic systems and 
suggests that rather than looking for the unpredictable 
behavior of individual trajectories, it is better to look for 
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more reliable and stable properties derived by averaging 
over ensembles of initial conditions. Prominent among 
such properties is the distribution of lifetimes, obtained 
from many runs with similar but not identical initial con- 
ditions. The theoretical prediction for a hyperbolic sad- 
dle is that the probability of decay is constant in time and 
independent of when the flow was started, giving for the 
distribution of lifetimes an exponential, as in radioactive 
decay [13, [H, fl^. Other functional forms are possible 
as well (see e.g.. |2Cll.[2lj|). but for the most part observa- 
tions in transitional shear flows are compatible with an 
exponential Q, [n], [H, [H, H, HI] . 

An exponential distribution is characterized by a char- 
acteristic decay rate or a characteristic lifetime t which 
is the time interval over which the survival probability 
drops by 1/e. How this lifetime varies with Reynolds 
number is currently under debate [2^ . [2^ [27| . If r di- 
verges at a finite Reynolds number, there is a critical 
value RCc above which turbulent flow does not relaminar- 
ize but persists forever. Such a divergence would imply 
that the system undergoes a transition from a transient 
chaotic saddle to a permanently living chaotic attractor 
in some form of 'inverse boundary crisis' [28j . However, 
if T does not diverge, turbulence in a pipe remains tran- 
sient for all Re. The chaotic saddle does not close to form 
an attractor and the turbulent 'state' stays dynamically 
connected to the laminar profile even at Reynolds num- 
bers higher than the ones where 'natural transitions' are 
reported to occur. This might open up new avenues for 
controlling turbulent motion. 

The prediction of an exponential distribution of life- 
times is an asymptotic one, valid for long times. On 
short times the distributions may follow a different func- 
tional form, as evidenced by the non-exponential parts 
in almost all distributions published so far. Moreover, 
the results may depend on additional parameters, such 
as an aspect ratio or the length of the pipe. The depen- 
dence on these parameters has not been studied so far. 
It is our purpose here to discuss some of these effects for 
transitional pipe flow. 
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We begin in section |TT] with a survey of previous experi- 
mental and numerical results. Section lllll then is devoted 
to an analysis of three effects: the dependence on the en- 
semble of initial conditions in section UlI Bl the variation 
of the characteristic lifetime with Re in IIII CI and the 
variation with the length of the pipe in section llll Dl We 
conclude with a summary and outlook in section IIVI 



II. SURVEY OF RESULTS 

As usual, the mean downstream velocity (m}, the di- 
ameter D of the pipe and the viscosity v of the fluid can 
be combined into the dimensionless Reynolds number, 

Re=M^. (1) 

V 

The pipe diameter D and the velocity (u) then define a 
unit of time D/ {u). Since the flow moves downstream 
with the mean velocity (w), time can be translated into 
distance traveled, so that the distance in units of the 
diameter equals the time t in units oi D/ (u) . Because 
of this relation between length and time, it is crucial to 
work with very long pipes so that the observation times 
become as large as possible. 

When the flow becomes turbulent, the friction factor 
increases. Therefore, either the forcing (pressure drop) 
has to increase so as to maintain the mean flow speed, 
or the mean flow speed will decrease, perhaps reducing 
the Reynolds number so much that the flow relaminar- 
izes [2^. Thus, many modern experiments work with a 
constant flow rate [lO], or with very long pipes [l^], in 
which the change in Reynolds number becomes negligi- 
ble as long as the turbulence remains confined to a small 
section of the pipe: In the range of Reynolds numbers 
studied here, the turbulence is localized in a region of 
about SOD length To measure lifetime statistics one 
can either follow a puff on its journey down the pipe and 
determine the downstream position where it decays; or 
one can choose a fixed downstream position, which corre- 
sponds to a fixed lifetime, and measures the probability 
that puffs survive the journey down the pipe up to this 
chosen point. 

The first approach was chosen by MuUin in a recent 
series of experiments [H, [23, [H, [Slj inspired by the nu- 
merical studies in . In a first group of experiments [Slj 
the decay of the perturbation could be detected with a 
camera that traveled with the perturbation downstream. 
The length of the pipe allowed for a maximal observation 
time of 500 units. The flow was perturbed by injecting 
six jets of different amplitudes, and 40 to 100 indepen- 
dent repetitions were taken for each Reynolds number. 
The asymptotic regime of the distributions of lifetimes 
was found to follow a law 



with a characteristic lifetime r(Re) depending on Re and 
an initial offset io ~ 100 . . . 150 before which no de- 
cay was observed. The strong increase of the charac- 
teristic lifetime with Re lead to the conclusion that it 
diverges at a finite critical Reynolds number. For the 
critical Reynolds number they give in [3l| the values 
Rec = 1710 ± 10 and Re^ = 1830 ± 10, and in f23| the 
values Rec = 1695 ± 20 and Re^ = 1820 ± 20, for two dif- 
ferent kinds of perturbations, described as 'strong' and 
'weak' types of perturbation, respectively. In an effort to 
address the dependence on the type of initial perturba- 
tions, they performed a second experiment with a slightly 
different perturbation protocol |i26|: In order to obtain 
more generic initial conditions the system was started at 
a higher flow speed, a perturbation that triggered tur- 
bulence was introduced and then the Reynolds number 
was reduced to the one for which lifetime statistics were 
collected. This gives another sample of initial conditions 
but limits the remaining observation time to less than 
450. With such a perturbation the characteristic life- 
times were compatible with 

r(Re) cx (Re^ - Re)-i±° 02 , (3) 

but now with a different critical Reynolds number of 
Rcc = 1750 ± 10. 

It is difficult to model the perturbations induced by 
jets in numerical simulations f[2^). but it is relatively 
straightforward though time-consuming to imitate the 
second protocol of MuUin, where the initial conditions 
are taken from a turbulent flow at higher Reynolds num- 
bers. Willis and Kerswell [27] did just that for five differ- 
ent Reynolds numbers and concluded that Rcc = 1870, 
as suggested by some experiments. However, when the 
analysis of their data points is corrected as suggested in 
(32i] , the demonstration of a divergence is less convincing 
and the data become compatible with the results of Hof 
et al [2^. 

The experiments by Hof et al. [22] just mentioned use 
a different approach. In a pressure driven flow through 
a thin pipe of only 4 mm diameter but 30 m length they 
realized dimensionless observation times of up to 7500 
units. Since the flow could not be visualized, the time 
and position of decay could not be determined directly. 
However, a laminar and a turbulent patch in the flow 
can easily be distinguished once they leave the pipe, so 
that it is relatively easy and straightforward to deter- 
mine whether the flow has stayed turbulent until it exits 
the pipe. Therefore, they could determine the proba- 
bility to be turbulent after a time period given by the 
distance between the perturbation and the outlet, as a 
function of flow rate. This gives P{t, Re) as a function of 
Re for t fixed. Collecting data for different t then gives 
the parameters in the lifetime distribution including the 
Reynolds number dependence. For short times, the data 
are within the error bars of [2^ , but for longer times they 
deviate from the divergent behavior implied by In- 
stead, it was found that the lifetimes are well represented 
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by an exponential variation, 

T^^ (Re) oc exp(a + 6Re) (4) 
with a = 55.3 and b = -0.032. 



III. LIFETIME DISTRIBUTIONS AND THEIR 
PROPERTIES 

In this section we study Hfetime distributions in pipe 
flow within direct numerical simulations. Since the cal- 
culations are extremely time consuming, we will not aim 
to repeat the pufF simulations of (27i] , but rather focus on 
short, periodically continued pipe sections, and then dis- 
cuss how these results scale up to turbulence in regions of 
the length of turbulent puffs. In the next subsection we 
first discuss features of individual turbulent trajectories, 
before turning to the ensemble dependence of lifetime 
distributions, the variations with Reynolds number and 
the length dependence. 

Individual trajectories were generated using the pseu- 
dospectral DNS code develo ped in [3^ and already used 
in our previous studies [^■[l5l. l22l . [33 |. Simulations of elon- 
gated puffs with the determination of their travel velocity, 
envelope and internal dynamics are given in [3^ . The 
code uses Fourier modes in downstream and azimuthal 
direction and Chebyshev polynomials in the radial direc- 
tion, and a projection method to eliminate the pressure. 
The simulations on pipe segments presented in this sec- 
tion are carried out with n Fourier modes in azimuthal 
and m Fourier modes in downstream direction, where 

\n\/Njnax + \m\/Zmax < 1 with Nmax = 16 and Z^ax 

increasing from 14 for the 'short' pipe of length 5D to 
Zmax = 25 and Z^ax = 45 for the 'medium' (L = 91?) 
and 'long' {L = 15D) pipes, respectively. Consequently, 
we consider up to 33 Fourier modes in azimuthal direc- 
tion. In downstream direction up to 31 modes are con- 
sidered for the short, 51 for the medium and 91 for the 
longest pipe. We use 49 Chebyshev polynomials for the 
expansion in the radial direction. This moderate reso- 
lution results from a compromise of accurate represen- 
tation of the dynamics and maximum simulation speed, 
required for good statistics. 



A. Features of individual trajectories 

Consider a perturbation of the laminar Hagen- 
Poiseuillc flow applied at time t — 0. The evolution of the 
initial condition u{t — 0) can be followed in time until it 
decays or reaches the maximum integration time in a sim- 
ulation or leaves the pipe in an experimental setup. Fig.[T] 
shows the evolutions of four different but similar initial 
conditions. As an indicator for the turbulent intensity, 
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FIG. 1: (Color online) Simulated evolution of four initial dis- 
turbances at Re = 1900 in a periodic pipe of length L = 15D. 
Plotted is the kinetic energy of the deviation from the lam- 
inar profile that is stored in the streamwise varying Fourier 
modes, eq. [5] Initial conditions are constructed from a mod- 
ulated Zikanov mode discussed in the main text. The four 
initial conditions differ by less than 0.5% in energy content. 
The pertubations first grow in energy and show an overshoot 
before settling down to the turbulent state. They then sud- 
denly decay without any prior indication and the energy of 
the perturbation decays monotonically. The chosen criterion 
for decay is based on the energy threshold Esd < 5 x 10~^. 



we take the energy of the three-dimensional structures. 




where if^ ^ denotes the (n, rn)-Fourier mode if the per- 
turbation field V = u — 2(1 — r'^)ez is decomposed into 
Fourier modes in azimuthal (m) and axial (n) direction. 
The energy content of the streamwise modulated Fourier 
modes is normalized by the kinetic energy of the laminar 
profile. 

Since a flow fleld only asymptotically reaches the lam- 
inar proflle exactly, 'decay' is defined as reaching a sit- 
uation where perturbations of the laminar profile are so 
small that the further evolution follows an exponential 
drop off. Any perturbation is therefore characterized by 
a lifetime that slightly depends on the chosen criterion 
to detect being close to the laminar profile. Technically, 
one introduces a cut-off threshold either on the kinetic 
energy stored in the deviation from the laminar profile 
/ 2(1— r^)e2|pdF, or on the kinetic energy ^ stored 
in streamwise invariant Fourier modes only. The thresh- 
old on these energies is chosen such that the further evo- 
lution can be described by the linearized equations, so 
that the system cannot return to the turbulent dynam- 
ics. The lifetime then is defined as the time it takes to 
reach this target region around the laminar profile. 
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FIG. 2: (Color online) Spatial structure of the types of ini- 
tial perturbations considered. Left panel: In plane velocity 
components of Zikanov's nearly optimal growing mode. Axial 
velocity components vanish. The translational symmetry of 
Zikanov's mode is broken by applying a 2-dependent twist. 
Right panel: Cross-section of a velocity field from a turbu- 
lent run at Re = 2150. The vector plot indicates the in-plane 
motion. Color coding is used for the axial velocity relative to 
the laminar parabolic profile. 



B. Ensemble dependence 

Now consider an ensemble of several different but sim- 
ilar perturbations. The collection of individual lifetimes 
can be used to estimate the probability P{t) to still be 
turbulent at some time t. A chaotic saddle should give 
rise to exponential asymptotic tails of this distribution 
that are independent of the choice of initial conditions 
but characteristic for the saddle. 

One type of perturbations we consider here is a pair 
of vortices as in the optimally growing modes identified 
by Zikanov [sgI . [stI . [3a | . In order to break translational 
symmetry they are modulated in streamwise direction by 
applying a z-dependent twist: 

uo{r,(f,z) :=^uzik(^r ,ip + (fosin(^z^ , , (6) 

where uzik is Zikanov's mode and L is the length of the 
computational domain used in our direct numerical simu- 
lation. The spatial structure is presented in Fig. [21 Mim- 
icking experimental protocols, where the spatial struc- 
ture of the perturbation is prescribed by the setup, the 
ensemble of initial conditions is constructed by varying 
the amplitude of the twisted Zikanov mode. A second 
type of perturbation is a snapshot from a turbulent run 
at Re = 2150 which is scaled in energy i.e., in amplitude. 
A cross-section is shown in Fig. [51 

Fig. [31 shows the lifetime of the perturbation as a 
function of the initial energy Eq- Regions of small and 
smoothly varying lifetime are clearly separated from re- 
gions of longer fiuctuating lifetimes. In regions with 
short lifetimes the flow relaxes quickly to the laminar 
profile. Towards the boundaries of these regions the life- 
times increase quickly and reach plateaus at the maxi- 
mal integration time. Magnifications of the plateau re- 
gions show chaotic and unpredictable variations of life- 
times fTTJ . The cliff-structure in the lifetime is due to the 
geometric features of the basin of attraction of the 1am- 
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FIG. 3: (Color online) Lifetime of initial conditions at 
Re = 1900. The pipe is L = 15D long. The initial conditions 
are constructed from a turbulent snapshot of a Re = 2150 
run (Fig. [2]). It is scaled in kinetic energy Eq. The black 
crosses belong to an ensemble constructed by varying the en- 
ergy form Eo = 8x 10'^ to Eq = 8.4 x 10"^ in 100 equidistant 
steps. Additional ensembles with 100 samples each focus on 
regions of high probability to observe long turbulent tran- 
sients. The yellow triangles indicate an ensemble that ranges 
form Eo = 8.0 x 10"^ to Eq = 8.05 x 10"^ in 100 steps. The 
red squares are located in the range between Eq = 8.23 x 10"^ 
and Eo — 8.28 x 10"^. Finally the cyan circles reach from 
Eo = 8.28 X 10"^ to Eo = 8.38 x 10"^ The corresponding 
lifetime distributions for all ensembles are shown in Fig. U 



inar profile and has been seen in pipe flow 13411 and low- 
dimensional representations of shear flows [l4i [T^ [39|. 

A first ensemble of initial conditions is constructed 
by varying the energy from Eq — 8.0 x 10^^ to Eq — 
8.4 X 10"^ in 100 equidistant steps. Other ensembles 
are chosen such that they provide a higher resolution in 
initial energies for the regions where high lifetimes are 
expected. 

Fig. [H shows P{t) for fixed Re calculated for the dif- 
ferent ensembles of initial conditions. Each individual 
distribution is characterized by an initial offset to « 75 
where no trajectory decays. Then one observes a mid- 
dle part that both in length and functional form differs 
among the various ensembles before asymptotically ex- 
ponential tails are reached. Within the limits of statis- 
tical uncertainties the decay rates (i.e., the slopes in a 
semi-logarithmic plot) are independent of the ensemble 
of initial conditions and encode a characteristic feature 
of the chaotic saddle. 

In view of the cliff-structure in lifetime the origin of 
the non-universal middle part becomes obvious: Only 
trajectories starting from the regions of fluctuating life- 
time reach the chaotic saddle from which they can decay 
at a characteristic rate. Initial conditions from the re- 
gions in between directly decay without having reached 
the saddle which leads to the non-universal initial parts 
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FIG. 4: (Color online) Lifetime distribution for fixed param- 
eters (Re — 1900 and periodic domain of length L = 15D) 
for different ensembles of initial conditions. Presented is 
data obtained from an ensemble of 200 samples constructed 
from twisted Zikanov modes (green down-facing triangles). 
The other ensemble result from scaling a turbulent flow field 
from a Re = 2150 run. The black crosses indicate an en- 
semble constructed from scaling the energy of the field from 
So = 8 X 10"^ to Eo = 8.4 X 10"^ in 100 equidistant steps. 
The lifetime as a function of this energy shows typical folds 
(see Fig.|3]). Based on this fold structure remaining ensembles 
indicated by the yellow triangles, red squares and cyan circles 
are tailored towards reaching the saddle i.e., additional runs 
were chosen in the 'right' energy region. The specific ensem- 
bles containing 100 samples each are described in Fig. |3]where 
the same symbols are used. While the initial part of the dis- 
tributions depends both in length and functional form on the 
chosen ensemble, the slope of exponential tails is universal 
within the limits of statistical accuracy. 



of P{t) in Fig.m The different plots in Fig. |4] correspond 
to different ensembles chosen from all possible initial con- 
ditions presented in Fig. [3l 

We now focus on the initial offset time Iq. It corre- 
sponds to the smallest time it takes for an initial con- 
dition from the chosen ensemble to decay. It evidently 
depends on the initial perturbation. In particular, this 
time scale can in principle be arbitrarily close to zero 
by including an arbitrarily small perturbation to the lin- 
early stable laminar profile in the ensemble. However, 
'typical' perturbations used both in simulations and lab 
experiments are characterized by a typical initial forma- 
tion time to (with w 100 . . . 150 in experiments). This 
can be rationalized as follows: 

We first note that to is large compared to the Lyapunov 
time measured in the turbulent motion [llj | , which gives 
a typical timescale for the dynamical separation of neigh- 
boring trajectories on the saddle. However, a trajectory 
does not necessarily start on the chaotic saddle. Its ini- 
tial condition is a flow field that hopefully initiates tur- 
bulence i.e., the trajectory approaches the chaotic saddle. 



but typically it does not belong to the saddle itself. In 
addition a trajectory starting its decay from the saddle 
has to follow the evolution through state space until it 
ends up in the vicinity of the laminar profile. Conse- 
quently the offset to contains two parts: The formation 
time tfor required to reach the turbulent 'state' and the 
decay time tdcc it takes to finally reach the neighborhood 
of the laminar state after the decay has been initiated. 
This can also be directly observed in Fig. [1] The pre- 
sented energy traces consist of three parts: the initial 
energy growth, the chaotic fiuctuations indicating turbu- 
lent dynamics on the saddle, and the decay towards the 
laminar state. 

As discussed in Section |ll] the decay time can be esti- 
mated from experimental observations that a puff decays 
while traveling about 50D downstream [l^]. This trans- 
lates into a time t^cc ~ 50. Theoretically it follows from 
the mechanism of decay: It starts first with a reduction 
of transverse modulations leading to a break out of the 
regeneration cycle supporting turbulent motion j40| and 
then shows a viscous damping of streamwise deviations 
from from the laminar profile. 

The formation time ifor can be estimated from the ex- 
perimental observation that is takes tfor about 40 for 
a perturbation (jet injection [l2|) to develop into an 
equilibrium puff. It is also observed in simulated short 
periodic pipes: Taking an initial condition that con- 
tains pairs of counter-rotating vortices in axial direction 
such as Zikanov's almost optimally growing mode, these 
small perturbations grow in energy by generating strong 
streaks in a so called lift-up process driven by the non- 
normality of the Navier-Stokes operator [4]|. The gen- 
erated streaks then become unstable, transverse vortices 
appear and the fiow turns turbulent. The formation time 
can therefore be estimated by the growing period of a 
Zikanov mode and originates from the non-normal char- 
acter of the Navier-Stokes operator. 

The discussed mechanism is also present in experi- 
ments where the fiow is perturbed by injecting fiuid jets 
perpendicular to the pipe These jets give rise to 

pairs of counter-rotating vortices that will draw energy 
from the base profile and grow by the same mechanism 
observed in the dynamics of a streamwise independent 
Zikanov modes. The experimental observation agrees 
with typical formation times of equilibrium puffs in sim- 
ulations in a periodic pipe of length L — 50D where a 
localized form of Zikanov's mode was used as an initial 
condition. 

We thus conclude that the specific choice of initial con- 
ditions affects the functional form of the lifetime distri- 
bution for small times. Both the initial time offset tg 
and the middle part of the lifetime distribution depend 
on the chosen ensemble of initial conditions. Universal 
properties of the turbulent state are only encoded in the 
asymptotic tails, which result from trajectories that actu- 
ally reach the turbulent state before their decay. The ex- 
ponential form of these tails is compatible with a chaotic 
saddle in state space and characteristic decay rates can 
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be 'measured' by fitting exponentials to the asymptotic 
tails of the lifetime distributions. 

In order to probe the chaotic saddle, i.e., analyze the 
universal asymptotic part of the distribution and mini- 
mize the high computational costs at the same time, it is 
favorable to choose an ensemble of initial conditions that 
is likely to reach turbulence. Modulated Zikanov modes 
seem to be a good choice whereas independent snapshots 
from a turbulent run at slightly higher Re tend to de- 
cay directly. This can for example be observed by com- 
paring the lifetime distributions (Fig. |4]) of the ensem- 
bles constructed by scaling the turbulent flow field and 
the Zikanov mode respectively. Although snapshots from 
runs at higher Reynolds number are complex and appear 
to be 'turbulent' they need not be located close to the 
chaotic saddle in the state space of the system at slightly 
lower Re. Typical length and time scales of turbulent 
motion change with Re so that turbulent snapshots at 
one Reynolds number might not 'fit' the dynamics at an- 
other Re [131 • In contrast, in the case of Zikanov modes 
the flow has enough time to adapt to the dynamics of the 
Navier-Stokes equations. Moreover, the Zikanov mode 
shares a characteristic pair of streamwise vortices with 
the edge state of pipe flow [sl, S^l • Since the edge state 
is located in-between the laminar state and the chaotic 
saddle, an initial condition close to the edge state should 
be especially efficient in initiating turbulence. 

To summarize: The specific form of a lifetime distri- 
bution does not only depend on the system parameters 
such as Reynolds number and — in the case of a sim- 
ulation — length of a periodic domain and resolution 
of the numerical representation that completely define 
the dynamical system. P(t) also depends on the ensem- 
ble of initial conditions. The large initial offset due to 
transient growth of initial perturbations and the initial 
drop generated by trajectories that decay directly with- 
out reaching the chaotic 'state' are not universal. In view 
of the cliff-structure in lifetime that shifts with Reynolds- 
number even choosing the same ensemble of initial condi- 
tions does not prevent complicated variations of the non- 
universal parts of P{t) with Re. Only the exponential 
tails in the asymptotic regime of the lifetime distribu- 
tion carry information about the chaotic saddle and its 
characteristic decay rate. Consequently, long observa- 
tion times reaching into the asymptotic range and initial 
conditions that have a high probability to reach turbu- 
lent dynamics are needed in any study of characteristic 
lifetimes. 



C. Reynolds number dependence T(Re) 

Lifetime experiments were performed in periodic do- 
mains of length L = 5D, L = 9D and L = 15D for 
various Reynolds numbers. The ensemble of initial con- 
ditions was constructed from twisted Zikanov modes of 
varying amplitude. For each Reynolds number at least 
100 independent trajectories were integrated up to a 
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FIG. 5: (Color online) Probability P{t) to still be turbulent 
after some time t for the 'short' L = 5D pipe. For each 
Reynolds number (given in the legend) 100 independent initial 
conditions have been integrated up to a maximum integration 
time of 750. Exponential fits to the tails of the distribution 
are indicated by straight lines in the semi-logarithmic repre- 
sentation. The measured slopes are shown in Fig. [6] 



maximum integration time of 3000i?/J7c; = 750D/(m) 
which is about twice the observation time available in the 
Manchester pipe and a tenth of the maximal observation 
time in the discussed experiments by Hof. Characteristic 
lifetimes were extracted from the slopes of exponentially 
varying tails of the measured lifetime distributions. 

Fig. O shows the probability P{t) to still be turbulent 
at some time i as a function of this time for the 'short' 
L — 5D pipe. The data points lie on straight lines in a 
semi- logarithmic plot, clearly indicating an exponential 
variation for large times. The slopes of the indicated 
fits which correspond to 1/r with r the characteristic 
timescale of the decay process are plotted in Fig. [6] as a 
function of the Reynolds number. 

The characteristic lifetime r increases rapidly with Re 
which, in previous studies [Hi [H, [sH , led to the conclu- 
sion that it diverges at a finite critical Reynolds number 
Kcr like 



r(Re) 



1 



Rcc — Re 



(7) 



In a linear plot of the inverse lifetime 1/r as a func- 
tion of Re this would correspond to a linear variation 
that crosses zero at the critical Reynolds number. Our 
data do not support this scaling. Indeed, 1/t approaches 
zero as we increase Re but there is no indication of a 
divergence. Instead the data is compatible with an expo- 
nential scaling which corresponds to a straight line 
in a semi-logarithmic representation. The values for the 
parameters a and b are listed in table HI 

Consequently, there is no evidence for a transition from 
a chaotic saddle to a permanent chaotic attractor, at least 
not close to a Reynolds number of order 2000, where 
critical values have been reported previously. 
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FIG. 6: (Color online) Inverse characteristic lifetimes for 
the 'short' pipe of length L = 5D, as extracted from the ex- 
ponential fits of P{t) (see Fig. [5} as a function of Reynolds 
number both on a semi-logarithmic (main) and a linear scale 
(inset). The error bars result from extreme fits still compati- 
ble with the data. Thus, they reflect both intrinsic statistical 
uncertainties of the exponential fits and errors resulting from 
estimating the beginning of the asymptotic regime. The data 
is compatible with an exponential variation of T(Re) but does 
not support the linear scaling proposed in the literature be- 
fore (cf. inset). Thus, the data does not support a divergence 
of T i.e., a zero crossing of 1/r at a finite Reynolds number 
close to 2000. 

Comparing with experimental results ^23\ , our numeri- 
cal studies confirm exponential lifetime distributions and 
that the characteristic lifetime does not diverge but grows 
exponentially with the Reynolds number. However, the 
parameters of the exponential scaling law do not match 
quantitatively. This shortcoming is addressed in the next 
section, where wc discuss how the length of the periodic 
domain used in our simulation affects the statistics. 



D. Extensitivity of r(Re) 

Having found the exponential scaling of characteris- 
tic lifetimes both in experimental works and in numeri- 
cal simulation, one can quantitatively compare both sys- 
tems. The main difference between both considered sys- 
tems is that in an experiment one observes a localized 
turbulent puff traveling through a very long pipe whereas 
in short simulated periodic pipes not the full spatiotem- 
poral structure but only the internal dynamics of a puff 
is captured. There is no coexistence of a turbulent re- 
gion and laminar flow and no dynamics of the fronts of 
a turbulent puff. If the periodic domain becomes long 
compared to all internal scales of a turbulent puff, in- 
cluding its overall size, features of the experiment should 
be quantitatively recovered. For shorter computational 
domains, however, finite-size effects are to be expected. 
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FIG. 7: (Color online) Lifetime distributions for the pipe of 
'medium' length L = 9D pipe. For each Reynolds number 
(legend) 100 trajectories have been analyzed. The straight 
lines indicate exponential fits to the tails of the distributions. 




t 



FIG. 8: (Color online) Same as Fig. [7] but for the 'long' pipe 
of length L — 15D. Note the long offset to especially for 
Re = 1900. 

We analyze the dependence on the length L of the 
computational domain by comparing the results from 
the short L = 5D pipe with additional simulations for 
a 'medium' L = 9D and a 'long' L — 15D pipe. The 
length L = 9D (rather than L = IQD) for the medium 
pipe was chosen such that periodic structures that might 
be favored by the periodicity of the small reference calcu- 
lation do not exactly 'fit' the new period in downstream 
direction. 

Fig. [7] shows lifetime distributions based on 100 indi- 
vidual runs at every Reynolds number for the medium 
pipe. Exponential fits to the tails are presented as 
straight lines in the semi- logarithmic plot. 

Fig.[S]shows the same data also based on the analysis of 
100 runs at each Reynolds number for the long periodic 
pipe. In this plot, the different parts of P{t) i.e., the 
initial decay followed by asymptotic tails that have to 
be used for measuring the characteristic lifetimes, are 
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FIG. 9: (Color online) Inverse characteristic lifetimes ex- 
tracted from simulations in periodic domains of different 
length (L = 5D (red squares), L = 9D (green circles), L — 
15D (blue triangles)) and from the experiment by Hof [22,] 
(black crosses) in a linear (left panel) and semi-logarithmic 
plot (right panel). The straight lines are exponential fits to 
the data, indicating that all data sets are individually com- 
patible with an exponential scaling of lifetime with Reynolds 
number. The slope of the exponential grows with the length 
of the periodic domain and approaches quantitatively the ex- 
perimental results. By extrapolation shown in Fig. [12] one 
can speculate that numerical and experimental results should 
match at L = 3073, which is about the length of an equilib- 
rium puff. 



obvious. 

In Fig. [5] the extracted inverse characteristic lifetimes 
from Fig. [7| and Fig. [5] are presented together with the 
L = 5D reference data. 

All three presented data sets show no evidence for a 
divergence. Each one is fully compatible with an expo- 
nential scaling of t as a function of Re. However the 
slope of the exponential is not universal but varies with 
the length of the computational domain. Large changes 
in the characteristic lifetime occur in a smaller interval of 
Reynolds numbers for a longer pipe. For Reynolds num- 
bers larger than about 1850 the characteristic lifetime 
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FIG. 10: (Color online) Autocorrelation functions for the 
downstream (Czz), azimuthal {C,pip) and radial (Crr) veloci- 
ties along the pipe axis. The correlations are evaluated in the 
comoving frame of reference of a turbulent puff [S^ and are 
normalized to one for vanishing axial shift Az. 

grows with the length of the computational domain. 

Hence, the characteristic lifetime is not a purely in- 
tensive measure. It scales with the size of the system 
under consideration. Such a scaling is compatible with 
the reasoning that has been developed for spatially ex- 
tended transient chaos, when the domains considered are 
larger than the correlation length. 

At the Reynolds numbers considered, the flow in the 
pipe is correlated in azimuthal, radial and also axial di- 
rection. The axial auto-correlation functions for the three 
velocity components are shown in Fig. [TOl They fall off 
rather quickly within about 2R. The short axial corre- 
lation length suggests that useful information about the 
interior dynamics of the flow can be obtained by studying 
relatively short domains. Specifically, for a length of lOi? 
the axial correlations in the downstream fluid are down 
to 30 per cent, but the computational advantages are 
enormous and allow for detailed studies of deterministic 
f43| and statistical properties [l5l|. 

A turbulent puff, on the other hand, defined via the 
total energy content, and averaged over time, is much 
longer, see Fig. [TlJ Therefore, it is plausible to think 
of the flow in a turbulent puff as succession of segments 
of length D, which are more or less uncorrelated. Let r 
be the probability for one of them to decay over a time 
interval 5t. Then the probability for N of them to decay 
simultaneously is . It has to be simultaneous, for oth- 
erwise a turbulent cell could main turbulence and trigger 
a spreading along the axis. Since r — 1/r, this implies 
that the lifetime rjv for the system of N cells is tjv = . 
With the representation ^ this implies a linear increase 
of the parameters a and b with the number of cells, i.e. 
the length of the turbulent section. This is nicely borne 
out by the data in Fig. [T^ 

Comparing with the characteristic lifetimes extracted 
from the experiments by Hof included in Fig. [9] the life- 
times extracted from simulations quantitatively approach 
experimental values as one increases the length of the 
computational domain. 
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TABLE L Parameters for the exponential law Q. The top 
three rows are from numerical simulations of periodic sections 
of length L = 5D, 9D and 15-D, respectively. The next to last 
line gives the value extrapolated for L = SOD and the last row 
gives the experimental data of p^ . 



FIG. 11: (Color online) Energy content of a turbulent puff, 
averaged over time in a co-moving frame (as discussed in [35l | ') 
in a logarithmic representation (main panel) and a linear on 
(inset). The rapid, exponential drop of underlines the strong 
localization properties of a turbulent puff and also shows that 
its length is rather well defined to be about 30 D. The spot 
moves from left to right. 
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FIG. 12: (Color online) The variation of inverse lifetime 
shown in Fig. |9] may be described as = exp(a + 6 Re) 
with two parameters a and 6. We plot these parameters as a 
function of the length L of the computational domain. The 
red boxes correspond to the parameter a/ 10 and the yellow 
circles to 1006. The experimental values from '22] are indi- 
cated by the green dashed horizontal lines. Linear fits to the 
data points (blue dashed lines) reach the experimental values 
at approximately the same length L « 3073 which is about 
the length of an equilibrium puff. 



The individual exponential variation of the lifetime 
with Re may be described using (jj]) with two parame- 
ters a and b which depend on the length L of the com- 
putational domain. In Fig. [12] we plot these parameters 
as a function of L. Linear extrapolation shows that both 
parameters independently approach the experimental val- 
ues - included as horizontal lines in Fig. [12]- at the same 
length. A quantitative match between experimental and 



numerical results can thus be expected at a length of 
L = 30D, which is approximately the length of equi- 
librium puffs observed in long pipes. Hence, the data 
suggests that the characteristic lifetime is an extensive 
quantity scaling with the length of a turbulent patch. 

Our findings support the assumption that turbulent 
puffs do not decay from their boundaries but decay is 
initiated in the bulk region. Short periodic pipes appear 
to faithfully reproduce the internal dynamics of a tur- 
bulent puff and the dynamics of the fronts is of minor 
importance. 



IV. CONCLUSIONS 

The results presented here confirm that lifetime dis- 
tributions of turbulent pipe flow asymptotically follow 
exponentials. While there may be differences for short 
times, the long time behavior is robust. This observation 
supports the idea that turbulent motion is generated by 
a chaotic saddle in state space. Features of the measured 
probability functions can be explained in terms of an en- 
semble of initial conditions that either directly decay or 
reach the strange chaotic saddle. From this saddle, tra- 
jectories decay at a constant escape rate, i.e., indepen- 
dent of the previous state. This situation is analogous 
to that of a particle moving in a com plex box with some 
holes through which it can decay [Til. l45j| , or an unstable 
nucleus subject to radioactive decay. After escaping from 
the saddle trajectories then follow the slow dynamics to- 
wards the linearly stable laminar profile. 

A thorough analysis of lifetime distributions from both 
experimental and numerical studies of pipe flow as well 
as plane Couette flow shows that the characteristic life- 
times grow with Reynolds number, but that they do not 
diverge at a finite value of Reynolds numbers. As a con- 
sequence, there is no evidence that the chaotic saddle in 
state space turns into an attractor by some sort of 'in- 
verse boundary crisis' [isj . Even for Reynolds numbers 
exceeding 2000 turbulent signals consist of transients that 
decay finally, though the time for decay can be very long. 
As regards much higher Reynolds numbers and in par- 
ticular the transition from puffs to slugs which increase 
in length, the present results suggest that also the decay 
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time T{Re) may become longer. But this is insufficient to 
suggest a transition to a permanent attractor, so that the 
question about a global bifurcation that turns the turbu- 
lent dynamics into an attractor remains open. In any 
case, the fact that turbulent motion stays dynamically 
connected with the laminar profile at Reynolds numbers 
exceeding 2000 suggests that turbulent fiow could be in- 
tentionally laminarized at minimal energetic costs. 

An intriguing finding of the present study is a depen- 
dence of the characteristic lifetimes on the spatial exten- 
sion of a turbulent region. The larger it is the less likely 
it decays. For the case of puffs, it is possible to extrap- 
olate the present numerical results to the experimental 
ones. The results for pipe flow are consistent with theo- 



retical models for spatially extended systems. It will be 
interesting to see to which extend the extensive scaling 
of lifetimes seen here in pipe flow is generic and can be 
found in other linearly stable flows such as plane Couette 
flow. 
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